\(\int (a+a \sec (e+f x))^m \, dx\) [345]

   Optimal result
   Rubi [A] (verified)
   Mathematica [B] (warning: unable to verify)
   Maple [F]
   Fricas [F]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 12, antiderivative size = 83 \[ \int (a+a \sec (e+f x))^m \, dx=\frac {\sqrt {2} \operatorname {AppellF1}\left (\frac {1}{2}+m,\frac {1}{2},1,\frac {3}{2}+m,\frac {1}{2} (1+\sec (e+f x)),1+\sec (e+f x)\right ) (a+a \sec (e+f x))^m \tan (e+f x)}{f (1+2 m) \sqrt {1-\sec (e+f x)}} \]

[Out]

AppellF1(1/2+m,1,1/2,3/2+m,1+sec(f*x+e),1/2+1/2*sec(f*x+e))*(a+a*sec(f*x+e))^m*2^(1/2)*tan(f*x+e)/f/(1+2*m)/(1
-sec(f*x+e))^(1/2)

Rubi [A] (verified)

Time = 0.07 (sec) , antiderivative size = 83, normalized size of antiderivative = 1.00, number of steps used = 3, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.250, Rules used = {3864, 3863, 141} \[ \int (a+a \sec (e+f x))^m \, dx=\frac {\sqrt {2} \tan (e+f x) (a \sec (e+f x)+a)^m \operatorname {AppellF1}\left (m+\frac {1}{2},\frac {1}{2},1,m+\frac {3}{2},\frac {1}{2} (\sec (e+f x)+1),\sec (e+f x)+1\right )}{f (2 m+1) \sqrt {1-\sec (e+f x)}} \]

[In]

Int[(a + a*Sec[e + f*x])^m,x]

[Out]

(Sqrt[2]*AppellF1[1/2 + m, 1/2, 1, 3/2 + m, (1 + Sec[e + f*x])/2, 1 + Sec[e + f*x]]*(a + a*Sec[e + f*x])^m*Tan
[e + f*x])/(f*(1 + 2*m)*Sqrt[1 - Sec[e + f*x]])

Rule 141

Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_))^(p_), x_Symbol] :> Simp[(b*e - a*f
)^p*((a + b*x)^(m + 1)/(b^(p + 1)*(m + 1)*(b/(b*c - a*d))^n))*AppellF1[m + 1, -n, -p, m + 2, (-d)*((a + b*x)/(
b*c - a*d)), (-f)*((a + b*x)/(b*e - a*f))], x] /; FreeQ[{a, b, c, d, e, f, m, n}, x] &&  !IntegerQ[m] &&  !Int
egerQ[n] && IntegerQ[p] && GtQ[b/(b*c - a*d), 0] &&  !(GtQ[d/(d*a - c*b), 0] && SimplerQ[c + d*x, a + b*x])

Rule 3863

Int[(csc[(c_.) + (d_.)*(x_)]*(b_.) + (a_))^(n_), x_Symbol] :> Dist[a^n*(Cot[c + d*x]/(d*Sqrt[1 + Csc[c + d*x]]
*Sqrt[1 - Csc[c + d*x]])), Subst[Int[(1 + b*(x/a))^(n - 1/2)/(x*Sqrt[1 - b*(x/a)]), x], x, Csc[c + d*x]], x] /
; FreeQ[{a, b, c, d, n}, x] && EqQ[a^2 - b^2, 0] &&  !IntegerQ[2*n] && GtQ[a, 0]

Rule 3864

Int[(csc[(c_.) + (d_.)*(x_)]*(b_.) + (a_))^(n_), x_Symbol] :> Dist[a^IntPart[n]*((a + b*Csc[c + d*x])^FracPart
[n]/(1 + (b/a)*Csc[c + d*x])^FracPart[n]), Int[(1 + (b/a)*Csc[c + d*x])^n, x], x] /; FreeQ[{a, b, c, d, n}, x]
 && EqQ[a^2 - b^2, 0] &&  !IntegerQ[2*n] &&  !GtQ[a, 0]

Rubi steps \begin{align*} \text {integral}& = \left ((1+\sec (e+f x))^{-m} (a+a \sec (e+f x))^m\right ) \int (1+\sec (e+f x))^m \, dx \\ & = -\frac {\left ((1+\sec (e+f x))^{-\frac {1}{2}-m} (a+a \sec (e+f x))^m \tan (e+f x)\right ) \text {Subst}\left (\int \frac {(1+x)^{-\frac {1}{2}+m}}{\sqrt {1-x} x} \, dx,x,\sec (e+f x)\right )}{f \sqrt {1-\sec (e+f x)}} \\ & = \frac {\sqrt {2} \operatorname {AppellF1}\left (\frac {1}{2}+m,\frac {1}{2},1,\frac {3}{2}+m,\frac {1}{2} (1+\sec (e+f x)),1+\sec (e+f x)\right ) (a+a \sec (e+f x))^m \tan (e+f x)}{f (1+2 m) \sqrt {1-\sec (e+f x)}} \\ \end{align*}

Mathematica [B] (warning: unable to verify)

Leaf count is larger than twice the leaf count of optimal. \(711\) vs. \(2(83)=166\).

Time = 7.36 (sec) , antiderivative size = 711, normalized size of antiderivative = 8.57 \[ \int (a+a \sec (e+f x))^m \, dx=\frac {30 \operatorname {AppellF1}\left (\frac {1}{2},m,1,\frac {3}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right ) \cos ^2\left (\frac {1}{2} (e+f x)\right ) \cos (e+f x) (a (1+\sec (e+f x)))^m \sin (e+f x) \left (3 \operatorname {AppellF1}\left (\frac {1}{2},m,1,\frac {3}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right )-2 \left (\operatorname {AppellF1}\left (\frac {3}{2},m,2,\frac {5}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right )-m \operatorname {AppellF1}\left (\frac {3}{2},1+m,1,\frac {5}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right )\right ) \tan ^2\left (\frac {1}{2} (e+f x)\right )\right )}{f \left (45 \operatorname {AppellF1}\left (\frac {1}{2},m,1,\frac {3}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right )^2 \cos ^2\left (\frac {1}{2} (e+f x)\right ) (1+2 m-2 m \cos (e+f x)+\cos (2 (e+f x)))+6 \operatorname {AppellF1}\left (\frac {1}{2},m,1,\frac {3}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right ) \sin ^2\left (\frac {1}{2} (e+f x)\right ) \left (-5 \operatorname {AppellF1}\left (\frac {3}{2},m,2,\frac {5}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right ) (1+2 m-2 (2+m) \cos (e+f x)+\cos (2 (e+f x)))+5 m \operatorname {AppellF1}\left (\frac {3}{2},1+m,1,\frac {5}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right ) (1+2 m-2 (2+m) \cos (e+f x)+\cos (2 (e+f x)))-48 \left (2 \operatorname {AppellF1}\left (\frac {5}{2},m,3,\frac {7}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right )-2 m \operatorname {AppellF1}\left (\frac {5}{2},1+m,2,\frac {7}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right )+m (1+m) \operatorname {AppellF1}\left (\frac {5}{2},2+m,1,\frac {7}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right )\right ) \cot (e+f x) \csc (e+f x) \sin ^4\left (\frac {1}{2} (e+f x)\right )\right )+40 \left (\operatorname {AppellF1}\left (\frac {3}{2},m,2,\frac {5}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right )-m \operatorname {AppellF1}\left (\frac {3}{2},1+m,1,\frac {5}{2},\tan ^2\left (\frac {1}{2} (e+f x)\right ),-\tan ^2\left (\frac {1}{2} (e+f x)\right )\right )\right )^2 \cos (e+f x) \sin ^2\left (\frac {1}{2} (e+f x)\right ) \tan ^2\left (\frac {1}{2} (e+f x)\right )\right )} \]

[In]

Integrate[(a + a*Sec[e + f*x])^m,x]

[Out]

(30*AppellF1[1/2, m, 1, 3/2, Tan[(e + f*x)/2]^2, -Tan[(e + f*x)/2]^2]*Cos[(e + f*x)/2]^2*Cos[e + f*x]*(a*(1 +
Sec[e + f*x]))^m*Sin[e + f*x]*(3*AppellF1[1/2, m, 1, 3/2, Tan[(e + f*x)/2]^2, -Tan[(e + f*x)/2]^2] - 2*(Appell
F1[3/2, m, 2, 5/2, Tan[(e + f*x)/2]^2, -Tan[(e + f*x)/2]^2] - m*AppellF1[3/2, 1 + m, 1, 5/2, Tan[(e + f*x)/2]^
2, -Tan[(e + f*x)/2]^2])*Tan[(e + f*x)/2]^2))/(f*(45*AppellF1[1/2, m, 1, 3/2, Tan[(e + f*x)/2]^2, -Tan[(e + f*
x)/2]^2]^2*Cos[(e + f*x)/2]^2*(1 + 2*m - 2*m*Cos[e + f*x] + Cos[2*(e + f*x)]) + 6*AppellF1[1/2, m, 1, 3/2, Tan
[(e + f*x)/2]^2, -Tan[(e + f*x)/2]^2]*Sin[(e + f*x)/2]^2*(-5*AppellF1[3/2, m, 2, 5/2, Tan[(e + f*x)/2]^2, -Tan
[(e + f*x)/2]^2]*(1 + 2*m - 2*(2 + m)*Cos[e + f*x] + Cos[2*(e + f*x)]) + 5*m*AppellF1[3/2, 1 + m, 1, 5/2, Tan[
(e + f*x)/2]^2, -Tan[(e + f*x)/2]^2]*(1 + 2*m - 2*(2 + m)*Cos[e + f*x] + Cos[2*(e + f*x)]) - 48*(2*AppellF1[5/
2, m, 3, 7/2, Tan[(e + f*x)/2]^2, -Tan[(e + f*x)/2]^2] - 2*m*AppellF1[5/2, 1 + m, 2, 7/2, Tan[(e + f*x)/2]^2,
-Tan[(e + f*x)/2]^2] + m*(1 + m)*AppellF1[5/2, 2 + m, 1, 7/2, Tan[(e + f*x)/2]^2, -Tan[(e + f*x)/2]^2])*Cot[e
+ f*x]*Csc[e + f*x]*Sin[(e + f*x)/2]^4) + 40*(AppellF1[3/2, m, 2, 5/2, Tan[(e + f*x)/2]^2, -Tan[(e + f*x)/2]^2
] - m*AppellF1[3/2, 1 + m, 1, 5/2, Tan[(e + f*x)/2]^2, -Tan[(e + f*x)/2]^2])^2*Cos[e + f*x]*Sin[(e + f*x)/2]^2
*Tan[(e + f*x)/2]^2))

Maple [F]

\[\int \left (a +a \sec \left (f x +e \right )\right )^{m}d x\]

[In]

int((a+a*sec(f*x+e))^m,x)

[Out]

int((a+a*sec(f*x+e))^m,x)

Fricas [F]

\[ \int (a+a \sec (e+f x))^m \, dx=\int { {\left (a \sec \left (f x + e\right ) + a\right )}^{m} \,d x } \]

[In]

integrate((a+a*sec(f*x+e))^m,x, algorithm="fricas")

[Out]

integral((a*sec(f*x + e) + a)^m, x)

Sympy [F]

\[ \int (a+a \sec (e+f x))^m \, dx=\int \left (a \sec {\left (e + f x \right )} + a\right )^{m}\, dx \]

[In]

integrate((a+a*sec(f*x+e))**m,x)

[Out]

Integral((a*sec(e + f*x) + a)**m, x)

Maxima [F]

\[ \int (a+a \sec (e+f x))^m \, dx=\int { {\left (a \sec \left (f x + e\right ) + a\right )}^{m} \,d x } \]

[In]

integrate((a+a*sec(f*x+e))^m,x, algorithm="maxima")

[Out]

integrate((a*sec(f*x + e) + a)^m, x)

Giac [F]

\[ \int (a+a \sec (e+f x))^m \, dx=\int { {\left (a \sec \left (f x + e\right ) + a\right )}^{m} \,d x } \]

[In]

integrate((a+a*sec(f*x+e))^m,x, algorithm="giac")

[Out]

integrate((a*sec(f*x + e) + a)^m, x)

Mupad [F(-1)]

Timed out. \[ \int (a+a \sec (e+f x))^m \, dx=\int {\left (a+\frac {a}{\cos \left (e+f\,x\right )}\right )}^m \,d x \]

[In]

int((a + a/cos(e + f*x))^m,x)

[Out]

int((a + a/cos(e + f*x))^m, x)